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Abstract 

An orbital current mode peculiar to deformed quantum dots is theoretically 
investigated; first by using a simple model that allows to interpret analyti- 
cally its main characteristics, and second, by numerically solving the micro- 
scopic equations of time evolution after an initial perturbation within the 
time-dependent local-spin-density approximation. Results for different defor- 
mations and sizes are shown. 
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The study of collective excitations in semiconductor quantum dots is currently attract- 
ing much interest. Recent experiments, using resonant Raman scattering in two-dimensional 
GaAs-AlGaAs quantum dots have probed both charge-density and spin-density col- 

lective excitations, as well as single-particle excitations. From the theoretical side, charge- 
density excitations have been investigated since several years ago with different approaches 
such as Hartree 0, Hartree-Fock and density-functional theory In addition, a gen- 
eral scheme to describe spin-density excitations of both longitudinal and transverse charac- 
ter was presented in Refs. |^,|^. All these theoretical calculations applied the well known 
Random-Phase Approximation in circularly symmetric dots, for which the angular momen- 
tum selection rules can be exploited. 

Experiments on deformed nanostructures are nowadays providing very interesting new 
pieces of information (for instance, in Ref. [0 on ellipsoidal deformations). This has 
prompted, in Refs. UsHTTj, the extension of theoretical approaches to the symmetry- 
unrestricted situation. Motivated by this exciting new direction of the quantum dot field 
we report in this communication on a novel class of collective excitations. It involves the 
generation of orbital currents and is peculiar to deformed quantum dots. We will show that 
the orbital current excitation (OCE) is strongly connected to the quadrupole charge-density 
excitation (QCDE) and that its most clear signature is in the magnetic dipole strength (Ml). 
The characteristics of this mode will be investigated first analytically, by using a simplified 
model, and second, by solving numerically the time- dependent Kohn-Sham equations cor- 
responding to this particular motion. It is worth to point out that current modes similar 
to this one have been measured [|l^ and theoretically predicted [1^ in atomic nuclei, and 



are also expected to exist in deformed metal clusters and in the condensate of trapped 
bosons. In fact, while writing the present communication we have become aware of a mi- 
croscopic calculation for metal clusters, using an schematic Random-Phase Approximation, 
that has been published in Ref. |jl5| and one for Bose gases which is in preparation |jT6l . 



Let us assume a perfectly elliptic quantum dot, whose electron p and kinetic energy r 
densities are functions of the ellipse contour lines 

where and Ry are the elhpse radii. If R^ = Ry we obviously recover the circularly 
symmetric densities. A time- dependent displacement may be represented as a(t)u(r ), where 
a{t) is a time dependent parameter and u(r) gives the vector displacement at point r. The 
corresponding displacement operator that acting on the iV-electron ground state gives 
the displaced state \d) = I^|0), is 



V = exp 



i ^ 

—a{t)Y,{u{ri)-pi + H.a. ) 



(2) 



where pj is the i-th electron momentum operator and H.a. stands for the Hermitian adjoint 
operator. 

The OCE is represented by the following displacement field 
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u(r) = X r + r]V{xy) = {-y{l - t]) , x{l + rf)) 



(3) 



It is a combination of a rigid rotation with z axis and a quadrupole distortion, weighted 
with a parameter rj. This is a divergency free field V ■ u = and, by choosing 

^ Ry-R^ 

it can be shown that the density variations up to second order in a are 

(5p = 

2 

This result implies that the collective motion does not modify the electron density and 
therefore, it is not affected by the Coulomb interaction. In fact the cost in the Coulomb 
energy is minimized by adding the quadrupole term 'W{xy) in the displacement field @. 
The collective motion associated with this field modifies the kinetic energy density and thus 
it is an example of an elastic (or shear) mode exhibited by a Fermi system, that have been 
extensively studied in atomic nuclei . The electronic motion is in fact a collective rotating 
flow along the ellipse contour lines. 

The frequency of the OCE may be estimated assuming the oscillator formula cuqce = 



A 

M' 



where the restoring force k = 2E{u) is fixed by the energy change E{u) associated 
with the displacement field M = m J u^pdr is the collective mass parameter and m is 
the single-electron effective mass. One gets 
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M = Nm{r^){l-r]'^) , (6) 

where E^l is the unperturbed kinetic energy. If we further assume (r^) = ^r'^N and = 
^EpN (where ep and are the Fermi energy and the Wigner-Seitz radius, respectively) we 
finally get 



^OCE 




3 mr^ y^l — rj^ 



'L=N-'/' . (7) 



This simple expression tells us how the electronic OCE scales with iV, deformation of the 
ellipse 7], and electronic density (r^ parameter). 

Repeating a similar treatment with a pure quadrupole displacement u = V{xy) the 
frequency of the QCDE can be estimated 



l^QCDE 



v^^o , (8) 



where ujq is the average parameter of the external confining potentials in x and y directions, 
assumed of parabohc type, i.e., ojq = {u^ + t^j/)/2 (see below). 

From a physical point of view we expect both modes, OCE and QCDE, to manifest 
themselves in the response to the magnetic orbital dipole (Ml) operator fisLz, and to the 
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electric quadrupole (E2) operator xy. However, to really ascertain what is the relative weight 
of each mode in these channels we need to perform a more microscopic calculation. This is 
our goal in what follows. 

We will describe the time evolution within density functional theory, in the local spin- 
density approximation. The time dependent Kohn-Sham equations are solved by discretizing 
the xy plane in a uniform grid of equally spaced points and using the Crank- Nicholson ap- 
proximation. Of course, the unperturbed state is the Kohn-Sham ground state, numerically 
obtained by solving the static Kohn-Sham equations in the same grid by a steepest descent 
method. Technical details of the method can be found in \^JJ\. 



The set of single-particle orbitals {(y9j(r)} evolves in time as 

d 

'^Ol'^iriir, t) = h^[p, m] ipinir, t) , (9) 

where the spin index is rj =t, i, and total density and magnetization are given in terms of 
the spin densities Pr)(r) = J2i lv^ir?(r)P, P = + Pi and m = — pi, respectively. The 
Hamiltonian /i^ in equation contains, besides the kinetic energy, the confining potential 
y(con/)|-j,'j^ the Hartree potential f^^^(r) = / dr'p{r')/\r — r'\, and the exchange-correlation 
piece v^^^\r) = ■^Sxcip.rn). The exchange- correlation energy density 8xc{p,rn) has been 
described as in Refs. pHTOf. 

To model the elliptic quantum dots we will follow the prescription of Refs. and 
consider the confinement produced by anisotropic parabolas with parameters ujx and uoy in x 
and y directions, respectively. We define the ratio (3 = Uy/ux and fix the centroid with the 
Wigner-Seitz radius as lUq = In terms of these parameters, the external confining 

potential reads 

v'-'\-) = l-^l^^,i^' + PV). (10) 

As discussed in Ref. [0, this is a reasonable approximation to the real confining potential 
in vertical quantum dots with rectangular (mesa) structure. 

An initial perturbation, modelling the interaction with the physical probe, is needed in 
order to excite the system and monitor its time evolution. This is achieved by modifying 
the orbitals with the displacement operator (^. With our previous discussion, three natural 
options for the displacement field u(r) come inmediately to mind: a pure rotation (twist), a 
pure quadrupole distortion and the combination given in Eq. (^) (orbital distortion). In the 
latter case, we fix i] with the elliptic potential parameters by assuming that Rx/Ry = (3, i.e., 
the ratio between x and y radii is given by the inverse ratio of the corresponding parabola 
coefficients. By looking at the density contour lines we have checked that this assumption 
is well satisfied. It yields 

Figure 1 shows the results for the = 20 electron dot confined with = LSIoq and 
P = 0.75 (i.e., u;^=0.29 H*, u;y=0.22 H*). This figure nicely confirms the results anticipated 
with the analytical model. The orbital Ml strength is divided into two clear regions. One at 
high energy, which is associated with the QCDE, and one at low energy associated with the 
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OCE. The relative weigth given to both states is sensitively controlled with the parameter 
rj. When this is forced to zero (twist), the QCDE takes a large part of the Ml strength. 
But when t] is taken according to the system deformation, the OCE is the dominant mode. 
The decoupling is not perfect, as it is in the simple model, and the short period signal of 
the QCDE can be clearly seen, although with much lesser amplitude than for the twist 
excitation. The analytical expressions (|^) and (H) yield for this case the numerical values 

0. 066 H* and 0.36 H*, respectively, which agree quite reasonably with the microscopic result. 
The OCE is strongly fragmented because of Landau damping. This can be seen from Fig. 

1, since the OCE overlaps with peaks associated with single particle excitations, obtained 
by keeping /i^ fixed to its static value in Eq. (|^). By contrast, the QCDE is moved to higher 
energy by the residual interaction and is not fragmented. The lower panels of Fig. 1 show 
the Ml signal after an initial quadrupole distortion. In this case, the QCDE is the most 
dominant and we have used a logarithmic vertical scale to show that the OCE intensity is 
around 2% of the QCDE one. We conclude from Fig. 1 that the OCE produces a clear 
signature in the experimental orbital Ml strength of elliptical dots. 

We have peformed the same analysis with the xy signal, corresponding to the E2 channel, 
with similar conclusions, although the QCDE is more dominant in this case. The fragmen- 
tation patterns are the same as for the Ml results already discussed and the percentage of 
the OCE highest peak with respect to QCDE is ~ 1%, ~ 75% and ~ 1% for the twist, 
orbital and quadrupole distortions, respectively. 

The elastic behaviour of the OCE, in contrast to charge- density excitations, can be 
appreciated from the time evolution of the different contributions to the total energy after 
the initial perturbation. As seen from Fig. 2, after a rigid twist the system increases only its 
energy in the external field. This triggers the motion and, as a consequence of total energy 
conservation, the other energy contributions begin to increase at the expense of the external 
field term. With perturbation (|^) the situation is qualitatively different. At t = there is 
an important increase in both kinetic and internal Coulomb energies. The first is due to the 
elasticity of Fermi systems mentioned before, while the second appears because the quantal 
densities are not really uniform within the ellipsoid, but have shell oscillations. The rotation 
along the elliptic contour deforms this structure and thus increases the internal Coulomb 
energy. Again, because of energy conservation, once the motion has started all contributions 
fluctuate in a correlated manner. 

Figure 3 shows the evolution with deformation of the Ml strength, for = 20 and 
Tg = LSIoq, after an initial orbital displacement (^. In each case, the positions of the 
analytical energies and (H) are indicated by arrows. The energy dependence of the OCE 
is in qualitative agreement with Eq. ([^), decreasing from (3 = 0.5 to 0.875. However, formula 
(0) begins to deviate from the miscrocopic result for the highest deformation {(3 = 0.5). The 
QCDE does not change much with deformation since we keep the centroid ujq fixed. It can 
also be seen from this figure that when jS approaches unity, i.e., the dot becomes circular, 
the OCE carries less strength with respect to the QCDE. In fact, for j3 = 0.875 we have 
used logarithmic scale, since the OCE height is only ~1% of the QCDE one. 

In Fig. 4 we show the Ml strength for a fixed f3 and varying the number of electrons. 
There is a tendency to decrease the OCE energy with increasing size, correctly pointed by 
the analytical expression (^. It is also seen that the fragmentation of the OCE increases 
with size because of the increasing role of Landau damping. 
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To summarize, the existence of an intense, low lying OCE in the Ml strength of elliptical 
quantum dots has been pointed out. The importance of this mode with respect to the QCDE, 
as well as its evolution with dot deformation and size, have been discussed. Our analysis is 
based on the time evolution of the microscopic Kohn-Sham equations and on a simplified 
model that allows to obtain analytical solutions. The validity of the analytical expressions 
has been checked by comparison with the microscopic result. The microscopic calculation 
has also shown a sizeable fragmentation of this low energy mode, that we attribute to an 
important role of Landau damping. 

This work was supported in part by Grant No. PB95-0492 from CICYT, Spain. 
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FIGURES 



FIG. 1. Results for the time evolution of an elliptic dot with N = 20, /3 = 0.75 and = l.Slag, 
with the three different initial perturbations (rigid twist, orbital and quadrupole distortions). Left 
panels display the simulated orbital Ml signal as a function of time while right panels show the 
corresponding strength functions in arbitrary units. Middle right panel also shows the independent 
particle strength function (dashed line) and the position of the analytical approximations (0) and 
(^ with arrows. 

FIG. 2. Time dependence of the different energy variations after an initial twist and orbital 
perturbations: kinetic (solid), internal Coulomb (dash), external confining field (dash-dot), ex- 
change-correlation (dots) and total (dash-dot-dot). Notice that the sum of kinetic. Coulomb, 
external field and exchange-correlation contributions yield the total energy increment and remains 
constant in time. 

FIG. 3. Dependence of the OCE and QCDE with deformation (/3 parameter, see text) for 
= 20 electrons and = LSIcq. The arrows indicate the energies given by Eqs. (0) and (P). 

FIG. 4. Dependence with size of the OCE and QCDE strength for a fixed deformation of 
/5 = 0.75. As in Fig. 3, the arrows indicate the energies of the analytical expressions and (P). 
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